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1.  Introduction 

The  purpose  of  this  report  Is  to  outlltie  our  one- dimensional  transport  model  for 
cosmic  rays  propagating  in  the  Interstellar  medium  (ISH) .  The  equations 
described  here  have  been  used  to  prepare  a  numerical  model  of  galactic  cosmic - 
ray  transport  which  is  used  for  investigations  of  cosmic-ray  composition,  the 
sensitivity  of  cosmic-ray  abundances  to  cross  section  data  uncertainties,  and 
weak  reacceleration  in  the  interstellar  medium. 

The  basic  assumption  intrinsic  to  the  equations  described  here  is  that  cosmic 
rays,  most  of  which  seem  to  originate  in  our  galaxy,  while  traveling  through  the 
ISH  are  subject  to  a  number  of  proces.ses  which  tend  to  destroy  a  sense  of 
direction  of  propagation  with  regard  to  their  sources .  As  a  result ,  Instead  of 
dealing  with  three-dimensional  continuity  equations  which  imply  knowledge  of 
spatial  structure  and  source  distributions,  we  are  Invoking  a  one -dimensional 
"compressed"  view  of  cosmic -ray  diffusion.  This  view  is  exemplified  by  the 
treatment  of  "reflections"  at  the  "boundary"  of  the  galaxy  which  are  represented 
simply  as  a  probability  of  escape  from  the  galaxy,  or  escape  mean  free  path. 
Such  an  assumption  seems  to  adequately  represent  much  of  the  available  data  on 
cosmic-ray  composition,  but  may  be  of  limited  usefulness  In  representing  other 
features  of  cosmic  radiation  such  as  anisotropy. 

Within  the  one -dimensional  transport  model  there  is  a  further  assumption  that 
many  of  the  processes  necessary  to  interpret  cosmic-ray  composition  are 
insensitive  to  the  detailed  behavior  of  individual  particles.  Included  in  this 
assumption  are  the  use  of  interaction  mean  free  paths,  averaging  over  electron 
stripping  and  attachment  processes  for  the  estimate  of  electron  capture  decay 
rates,  and  representing  reacceleration  processes  with  a  Green's  function. 

The  processes  necessary  to  Interpret  cosmic-ray  data  are  approximated  using  a 
few  parameters  which,  It  is  hoped,  will  reveal  some  information  about  the 
structure  of  the  ISH  as  well  as  allowing  "background"  to  be  removed  in  estimating 
cosmic-ray  source  abundances.  The  parameters  of  principal  interest  are  the  mean 
amount  of  Interstellar  material  traversed,  the  mean  density  of  the  confinement 
region,  and  the  average  time  between  acceleration  at  the  source  and  loss  by 
escape  from  galactic  confinement.  These  parameters  are  well-known  and  often 
quoted  in  association  with  cosmic-ray  data  analysis. 

Additional  parameters  of  cosmic-ray  transport  are  discussed  iu  this  report. 
These  parameters  include  the  mean  strength  of  weak  shocks  in  the  ISH.  the  mean 
time  (or  pathlength)  between  encounters  with  these  shocks,  and  a  measure  of  the 
mean  inhoeogenelty  in  the  ISM.  These  parameters  introduce  additional  physical 
content  into  the  cosmic -ray  transport  equations  which  may  explain  certain 
anocaalles  in  the  cosmic-ray  composition  and  spectra.  The  additional  parameters 
may  also  be  helpful  in  ridding  the  cosmic-ray  transport  equations  of  ad  hoc 
parameters  and  functional  dependencies  which  have  been  Introduced  to  provide  good 
fits  to  the  cosmic -ray  data. 
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The  discussion  o£  this  report  is  arranged  into  sections.  Section  2  suniniarizes 
the  notation  used  in  this  report.  Each  variable  Is  identified  and  the  usual 
units  are  stated.  Section  3  contains  several  defining  and  kinematic  relations 
among  a  nusdier  of  the  variables.  Sections  4-10  each  describe  a  different 
constituent  part  of  the  transport  equation,  which  is  given  in  its  entirety  in 
Section  11.  Each  part  represents  the  infinitesimal  change  in  the  cosmic -ray  flux 
density,  Ji(E).  due  to  a  different  process.  Variovis  approximate  solutions  to  the 
cosmic-ray  transport  equations  are  possible  when  a  limited  number  of  constituent 
parts  are  combined  or  when  a  limited  number  of  cosmic-ray  species  are 
considered. ‘  Section  12  addresses  tdie  extension  of  the  transport  eqi;eations  to 
inhomogeneous  media. 


‘For  example,  G.  Gloecklcr  and  J.R.  Jo.ippi,  "Physical  Basis  of  the 
Transport  and  Composition  of  Cosmic  Rays  in  the  Galaxy,"  Phys.  Rev.  Lett.  22, 
1448-1453  (1969) 
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7.,  Hotatlon 


A* 

A, 


J(E) 

J(P) 

H 

m 

N(E) 

N(P) 

Na 

n 


n« 

P 


P 

Q(E) 

Q(P) 

q(E) 


Epioaw 


t 

V 


X 


isotopic  mass  number  of  a  cosmic-ray  species 

isotopic  mass  number  of  a  component  of  the  transport  medium 

speed  of  light  (cm  s'*) 

kinetic  energy  per  nucleon  (MeV) 

particle  flux  spectmm  (particles  cm'*  s'*  MeV'*) 

particle  flux  spectrum*  (particles  cm'*  s'*  [MeV/c]'*) 

mass  (MeV/c*) 

atonic  mass  unit  »  931.5  HeV/c* 

nxunber  density  spectrum  (particles  cm"*  MeV*) 

number  density  spectrum  (particles  cm**  [MeV/c]'*) 

Avogadro's  number 

number  density  of  transport  medium  (cm'*) 

number  density  of  a  component  of  the  transport  medium  (cm'*) 

momentum  (HeV/c) 

momentum  per  nucleon  (MeV/c) 

injection  rate  spectrum  (particles  cm"*  s'*  HeV*) 

injection  rate  spectmm  (particles  cm*  s'*  [MeV/c]'*) 

injection  spectrum  per  xmit  pathlength  (particles  cm'*  g*  MeV**) 

inverse  mean  free  time  or  rate  of  a  process  (s'*) 

inverse  mean  free  pathlength  of  a  process  (cm*  g'*) 

total  energy  (HeV) 

time  (independent  variable]  (s) 

velocity  (cm  s’*) 

pathlength  (independent  variable]  (g  cm'*) 
isotopic  charge  number  of  a  cosmic- ray  species 


y  «  shock  reacceleration  index 

p  «  saass  density  of  transport  medium  (g  cm'*) 

-  partial  mass  density  of  a  component  of  the  transport  medium  (g  cm"*) 
r  -  /naif life  for  nuclear  decay  (s) 

0  -  time  rate  of  energy  loss  per  nucleon  (HeV  s'*) 

w  -  stopping  power  per  nucleon  (MeV  cm*  g^) 


In  general,  toman  subscripts  label  cosmic-ray  species  and  greek  subscripts  label 
components  of  the  transport  medium  (XSH). 


*The  two  flux  spectra  mentioned  here  are  different  functions  and,  according 
to  coachematlcal  convention,  should  be  represented  as  such,  i.e.,  Jg(E)  and  Jp(p) . 
Since  we  will  have  no  occasion  to  mention  such  functional  fortes  as  Jj,(E) ,  we  have 
chosen,  for  the  sake  of  simplicity,  to  use  the  simplified  notation.  Such 
considerations  are  also  applicable  to  Che  number  density  spectrum.  M(E),  and  the 
injection  rate  spectrum,  Q(E). 
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3.  Rslstlons 

The  following  relations  exist  between  the  quantities  defined  in  Section  2 .  This 
is  not  intended  to  be  a  complete  or  alnimal  set  of  equations .  Several  equations 
serve  as  definitions  of  quantities  and  several  others  are  standard  results  of 
special  relativity;  the  remainder  may  be  derived  from  these. 


T*  -  +  mV 

13.11 

E  -  (T  -  Mc*)/A 

13.2] 

P  -  P/A 

13.3] 

m  ••  M/A 

13.4] 

dE  pc* 

dp  (B  *  me*) 

(3.5) 

X  -  pvt 

{3.6J 

d  3 

—  «•  pv  ■— 

dt  dx 

(3.7) 

J(E)  -  vN(E) 

(3.81 

N(£)  d£  -  M(p)  dp 

(3.9! 

J(E)  «  I5(p) 

(3.10] 

aNCE)  aj(E) 

at  dx 

(3.111 

Bp—  -  r>vr^ 

(3.12J 
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vQ<E)  -  QCP) 

P«  -  n»  A,  /  N* 

Loirents  time  dilation  factoir 

q(E)  -  Q(E)/pv 


(3.U1 

13.14] 

13.15] 

13.16] 
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4.  Source  Tern 

The  source  abundance  contribution  to  the  transport  of  the  1“  cosmic -ray  species 
Ifi  taken  to  be  of  the  form; 

aH,(p) 

- a(p)  [4.1] 

at 

where  Qj(p)  Is  the  rate  of  Injection  of  the  1*^  cosmic -ray  species  Into  the 
interstellar  medium  as  a  number  density  spectrum.  Often  It  is  assumed  that 
Q  a  where  T  &  2.* 

Using  [3.10]  and  [3.13]  the  contribution  to  the  1*^  component  of  the  cosmic-ray 
flux  spectrum  is: 

aJi(E) 

-  -  v<i(E)  [4.2] 

at 

Using  [3.7],  [4.2]  may  be  expressed  with  pathlength,  x,  as  the  independent 
variable : 

aJ,(E)  Q,(E) 

-  -  -  -  v<l(B)  [4.3] 

dx  p 

where  q,(E)  is  the  number  density  spectrum  of  cosmic  rays  injected  during  the 
period  that  the  I"*  species  of  cosmic  rays  of  energy,  B,  passes  tiicou^  a  unit 
pathlength  In  the  interstellar  medium  ([3.16]). 


*  Por  example,  J.F.  Ormes  and  R.J.  ProtJxeroe,  *lspIicatlons  of  H£AO  3  Bata 
for  the  Acceleration  and  Propagation  of  Galactic  Cosmic  Rays,”  Ap.  J.  272,  736- 
764  (1983) 
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5.  E«c«pe  Tern 

The  continuity  equation  nay  be  structured  to  approxiouite  diffusion  in  the  galaxy 
vlth  loss  at  the  boundaries  by  inclusion  of  a  galactic  escape  term.  A  single 
rate  of  escape  from  the  system  is  defined.  The  rate  may  be  dependent  on  the 
energy,  charge,  and  mass  of  the  transported  cosmic  ray.  This  transport  model 
Is  then  tensed  the  *leaky>box"  model. ^ 

The  one>diisenslonal  diffusion  equation  Is: 

aN(E) 

-  -  •R.CE)  N(E)  [5.1] 

dt 

where  R,(E)  is  the  time -rate  of  escape  from  the  galaxy,  and  time  and  pathlength 
are  related  by  (3.6).  For  convenience  we  have  dropped  the  species  index  from 
both  the  nu]id}er  density  spectrum  and  the  escape  rate;  it  is  understood  that  both 
quantities  may  depend  on  the  charge  and  mass  of  the  transported  species.  The 
species  Index  will  often  be  omitted  In  future  sections  when  no  change  in 
composition  occurs  for  a  particular  process. 

Applying  relations  [3.7],  [3.10],  and  [3.12]  as  in  Section  4  the  change  in  flux 
with  respect  to  pathlength  is: 

3J(E) 

- -  -  .r.(E)  J(E)  [5.2] 

dx 

where  r,(E)  is  the  inverse  mean  free  pathlength  for  galacclc  escape  as  a  function 
of  energy,  charge,  and  mass  of  the  transported  species. 

Kore  complex  diffusion  models  tsay  be  constructed  by  multiple  application  of 
{5.2}.  The  "nested  leaky-box"  is  an  example  of  one  such  model.  Additional 
transport  parameters  are  required  in  multiple  leaky  box  models. 

Pathlength  distributions  are  often  used  for  the  analysis  of  cosmic *ray  data.* 
Arbitrary  pathlength  distributions  are  not  generally  consistent  ^*lth  cosmic-ray 


'R.  Cowslk  aaid  L.W.  Wilson,  Proc.  13ch  Internat.  Cosmic  Ray  Conf.  I,  500 
<1973) 

*H.  Garcia-Kunos,  J.A.  Simpson,  T.G.  Guzik,  J.P.  Uetel,  and  S.H.  Margolis, 
“Cosmic-Ray  Propagation  in  the  Galaxy  and  in  th**  Heliosphere;  The  Path-Length 
Distribution  at  Lew  Exutv^,*  Ap.  J.  Suppl.  64,  ,^69-304  (1987) 
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diffusion  equations  when  energy  losses  and  gains  are  Inporc&nt.^  ’ih%re  is  no  use 
nade  of  pathlength  distributions  in  this  report.  To  achieve  the  effects  of  such 
distributions  on  cosaic>ray  composition,  one  must  introduce  additional  physical 
models  which  may  be  represented  as  terms  in  the  differential  equations  described 
in  this  report.  The  physical  models  are  likely  to  be  accompanied  by  free 
parameters  which  could  have  some  astrophysical  significance. 


*J.A.  Lesniak,  "The  Extension  of  the  Concept  of  the  Cosolc-Ray  Path-Length 
Distribution  to  Eonrelacivlstic  Energies,”  Ap.  Space  Sci.  63.  279-293  (1979) 
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6 .  Acceleration  Term 

A  term  In  the  transport  equation  for  acceleration  nay  be  derived  from  the 
assumptions  that  (1)  particles  are  conserved  In  the  acceleration  process,  (11) 
a  monochromatic  number  density  spectrum  Is  spread  into  a  power- law  spectrum  In 
momentiua,  p*^,  by  a  single  encounter  with  a  shock,  and  (111)  particles  are  not 
decelerated  by  this  process.^  Under  these  assumptions  the  following  result  Is 
derived: 


aN(p) 

at 


P 

-R.(P)  N(p)  +  (7-l)p^|  R,(p)  p''-‘  N(p)  dp 

0 


[6.1] 


where  R,  Is  the  inverse  mean  free  time  for  acceleration  events  (shock  encounter 
rate),  y  Is  the  shock  reacceleration  index,  and  p  Is  a  dummy  Integration  variable 
which  may  be  interpreted  as  the  cosmic-ray  momentum  prior  to  acceleration. 

In  terms  of  the  cosmic -ray  flux  spectrum  [6.1]  may  be  rewritten  using  [3.10]: 


aj(E) 

at 


-r.':e)  j(e) 


+ 


E 

I 

(7-l)p'^ 

« 

0 


R,(E) 


V‘  J(E)  d£ 


[6.2] 


where  p,  E,  and  v  are  dummy  quantities  with  the  same  relations  as  p,  E,  and  v. 
Converting  to  pathlength  as  Independent  variable  using  [3.7]  and  [3,12]: 


aJ(E) 


ax 


E 

-r,(E)  J(E)  +  (7-l)V‘p-^|  t,(E)  pT'-*  J(£)  d£ 

0 


[6.3] 


where  r,  is  the  inverse  mean  free  pathlength  for  acceleration  events.  The 
pcceleratlon  rates,  r,  and  R,,  can  be  species  dependent. 

By  addition  of  several  terms  of  the  form  of  the  right  hand  side  of  [6.3],  that 
aquation  may  be  generalized  to  Include  several  different  shock  strengths  each 
of  walch  Is  encountered  at  a  different  rate. 


’a.  Wandel,  D.  Elchler,  J.R.  Letaw,  R.  Silberberg,  and  C.H.  Tsao, 
“Distributed  Reacceleration  of  Cosmic  Rays,"  Ap.  J.  316,  676-690  (1987) 
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7.  lonlcatlon  U>ss  Tens 


As  cosmic  rays  move  through  the  ISM  some  of  their  energy  is  lost  due  to 
ionization  and  excitation  of  interstellar  matter.  In  general,  interstellar 
material  may  be  composed  of  several  atomic  species  labeled  by  k.  The  rate  of 
energy  loss  of  a  cosmic  ray  is  dependent  on  its  charge,  mass,  and  energy,  as  well 
as  the  composition  of  the  Interstellar  medium.  If  energy  losses  are  assumed  to 
occur  continuously,  the  ionization  loss  process  is  governed  by  the  following 
equation: 


3N,(E) 

at 


(P*/P)  Q*(E)  N,(E) 


[7.11 


where  Oj^fCE)  is  the  rate  of  energy  loss  per  nucleon  per  xmit  time  in  medium  k  for 
one  cosmic*ray  species.  The  factor  weights  tlie  stopping  powers  in  each 
component  of  the  ISM  according  to  their  abundance. 


In  terms  of  the  cosm'f'^-ray  flux  with  pathlength  as  an  independent  variable: 


aJ,(E) 

ax 


(P*/'P)  Ji(E) 


[7.2] 


where  «„^(E)  is  the  stopping  power  of  a  nuclide  in  energy  per  nucleon  per  unit 
pathlength  in  medium  x.  Heavy  ion  stopping  powers  have  been  reviewed  by  Ahlen.* 


'S.P.  Ahlen,  "Thecretical  and  experimental  aspects  of  the  energy  loss  of 
relativistic  heavily  ionizing  particles,"  Rev.  Hod.  Phys.  52,  121'173  (1980) 
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8.  Httclesr  Fxagnentation  Term 

The  change  in  nun&er  density ,  N,(E) ,  of  species  i  due  to  nuclear  fragmentation 
taking  place  in  a  medium  with  several  components  labeled  by  k  may  be  represented 
using  tlie  equation: 


at 


^  W*  n,(e) 

K 


00 


S<(E,E)^  Nj<£)  d£ 


[8.1] 


where  E  Is  the  energy  of  the  incident  cosmic  ray  (species  j)  and  E  is  the  energy 
of  the  resultiitg  spallation  fragment  (species  i) , 


Rt<E)«  -  lU  V  a(E)«,  [8.2] 

is  the  rata  of  destruction  of  species  1  in  medium  k  at  energy  E,  and  ^(E)^  is  the 
total  cross  section  for  species  change  (nearly  equal  to  the  total  Inelastic  cress 
section  at  higher  energies) ,  and 

Rf(E.£)4R  -  n*  V  a{E.E)^  (8.3) 

is  the  rate  of  production  of  species  i  at  energy  E  from  j  at  energy  iT  end 
velocity  v  in  medium  k,  and  o(E,E)^  is  the  partial  cross  section  for  production 
of  species  i  from  species  j .  Note  that  both  rates  for  nuclear  fragmentation  are 
explicitly  dependent  on  the  velocity  of  the  incident  nuclide  and  the  number 
density  of  the  ISH. 

In  terms  of  the  cosmic-ray  flux  spectrum  with  pathlength  as  the  independent 
variable  using  (3.8]  and  (3.11]: 


«0 


aj,(E) 

dx 


^  r,(E)«  J,(B) 


r,(E.E)^  Ji(£)  d£ 


[8.4] 


K  j  ,K  0 

where  the  r,  are  Inverse  mean  free  pathlengths  corresponding  to  the  rates,  R,, 
which  are  obtained  by  appl leaden  of  (3.12]  and  (3.14): 

r,(E)*  -  N,  ipjp)  o(E)«  /  A.  [3.51 


and 


r,(E.£)^  -  N*  Mp)  a(E.£)^  /  V 


[6.6] 
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where  is  Avogadro's  number,  is  the  partial  density  and  is  the  atomic  mass 
of  ISK  component  k.  We  note  that  pathlength  is  a  convenient  choice  of  the 
Independent  variable  in  cosmic -ray  transport  computations  because  the  r,  are  not 
explicitly  dependent  on  velocity  or  on  the  absolute  density  of  the  ISM. 

Tl\e  total  and  partial  cross  sections  for  heavy- ion  spallation  reactions  on 
components  of  the  interstellar  meditim,  especially  on  hydrogen  and  helium,  are 
an  Important  requirement  for  the  analysis  of  cosmic-ray  composition.  New 
measurements  of  spallation  cross  sections  are  reported  in  the  bi-annual 
Proceedings  rf  the  International  Cosmic  Ray  Conference.  Semi-empirical 
formulations  have  been  developed  by  Silberberg.  Tsao,  and  Letaw.’^'* 

Two  classes  of  cosmic-ray  heavy  ions  are  distinguished  according  to  their  history 
of  fragmentation  in  the  ISM.  Cosmic-ray  primaries  are  heavy  ions  which  have  not 
xragmented  in  the  ISM.  Cosmic -ray  secondaries  are  heavy  ions  which  have  arisen 
from  one  or  more  spallation  reactions  in  the  ISM.  Any  nuclide,  in  principle, 
can  arise  from  either  of  these  processes;  however,  H,  He,  C,  0,  Ne,  Mg,  Si,  and 
Fe  are  mostly  primary,  while  Li,  Be,  B,  F,  Na,  Al,  and  nuclei  with  21  £  Z  5  23 
are  mostly  secondary. 

The  ratio  of  stable  secondaries  to  primaries,  such  as  B/C,  is  determined  by 
oeas".rable  nuclear  fragmentation  cross  sections  and  the  escape  rate  from  the 
galaxy,  and  is  independent  o^  the  ISM  density  distribution.  Such  secondary  to 
primary  ratios  may  be  used  to  establish  aspects  of  the  Interstellar  diffusion 
process,  such  as  the  average  ament  of  matter  that  cosmic  rays  have  passed 
through. 


*R.  Silberberg  andC.H.  tsao.  ■Faiclal  Cross-Sections  In  High-Energy  Nuclec" 
Reactions,  and  Astrophyslcal  Applicaticns.  I.  Targets  with  Z  S  28.  li . 
Targets  Heavier  than  Nickel.,"  Ap.  J.  Suppl.  25,  313  and  333  (1973) 

*®J.R.  Letaw,  R.  Silberberg,  and  C.H.  Tsao,  "froton-Nucleus  Total  Inelastic 
Cross  Sections:  An  Empirical  Formula  for  £  >  10  KeV,”  Ap.  J.  Suppl.  SI,  271- 
276  (1983) 
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9.  Nuclemr  Decay  Term:  Particle  Eaiaalon 

The  following  equation  accounts  for  the  rate  of  change  of  a  cosmic>ray  nun&er 
density  spectrum  due  to  nuclear  decay  by  particle  emission: 


aN,(E) 

at 


-R,(E),  N.(E) 


Ro(E),  IS^(E) 


[9.1] 


where  S^CE),  is  the  rate  of  decay  of  species  i  at  energy  E  and  is  the  rate 

of  decay  of  species  j  into  species  i  at  energy  E. 


B^(E).  - 

In  2  /  (1  E/mc*)  r, 

[9.2] 

8^(E)^  - 

in  2  /  (1  +  E/mc*>  ry 

[9.3] 

where  r,  is  the  halflife  of  species  i  and  is  the  halflife  of  species  j  for 
decay  into  species  i.  All  rates  include  only  particle  emission  decay  modes; 
electron  capture  decay  modes  are  specifically  excluded  and  are  considered  in 
Section  10. 


In  terms  of  cosmic-ray  flux  with  pathlength  as  the  independent  variable  using 
[3.11]  and  [3.12]; 


5J,<E) 


dx 


-  -r,(E)t  J,(E)  +  )  r^(E),  J,(E) 


j 


[9.4] 


where  the  r„  are  inverse  mean  free  pathlengths  corresponding  to  the  R^. 

tjCE),  -  In  2  /  p  V  (1  +  E/ffic*)  r,  [9.5] 

and 

r^(E)^  -  In  2  /  p  V  (1  +  E/mc^)  .  [9.6] 


The  rate  of  decay  of  cosmlc-tay  species  depends  on  the  mean  density  of  the  ISK. 
If  the  mean  pathlength  in  the  ISH  has  been  established,  then  [3.6]  may  be  applied 
CO  estimate  the  mean  confinement  time.  Nuclides  such  as  “Al,  and  **C1  wltli 
half lives  of  about  one  million  years  have  been  analysed  for  this  purpose." 


"Fox  example,  H.E.  Wledenbeck,  "The  Isotopic  Composition  of  Cosmic  Ray 
Chlorine,"  Proc.  19th  Internet.  Cosmic  Ray  Conf.  2,  84-87  (1985) 
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10.  Huclear  Decay  Tera:  Electron  Capture 


Electron  capture  decay  of  cosqIc  rays  is  suppressed  because  hlgh>energy  heavy 
ions  passing  through  relatively  snail  anounts  of  material  are  likely  to  be 
stripped  of  their  electrons,  lover  energy  and  higher  charged  ioxis  may  have  a 
significant  probability  of  having  one  or  more  attached  electrons.  The  attached 
electrons  make  electron  capture  decay  possible  and  electron  capture  decay  may 
ultimately  have  an  observable  effect  on  composition.  The  importance  of  electron* 
capture  decay  modes  is  increased  in  distributed  acceleration  models  because  ions 
which  decay  freely  at  low  energies  are  mixed,  after  acceleration,  with  high* 
energy  components. 


The  nxinber  density  spectra  of  cosmic  rays  with  no  electrons  attached,  I^(E) ,  and 
with  one  electron  attached,  N‘(E},  are  governed  by  the  following  coiipled 
equations  if  the  probability  of  two  attached  electrons  is  negligible. ** 


aN«,(E: 

at 


R*(E)*  N®»(E)  .  B.(E)«  N\(E) 


R«(E)^  N‘j(E) 


tio.i] 


aN\(E) 

at 


^  R*(E)*  N«t(E)  *  B.(E)«  N‘i(E) 

K 


•  ME),  N\(E) 


(10.21 


where 

ME).  -  In  2  /  (1  +  E/mc»)  (O,  110.3] 

and 

ME),,  -  In  2  /  (1  +  E/mc»)  (10.4] 

are  the  rate  of  electron -capture  decay  of  species  I  at  energy  E,  and  the  rate 
of  electron-capture  decay  of  species  J  Into  species  i  at  energy  B,  respectively, 
and  the  electron  attachment  rate,  R.(E),  and  electron  stripping  rate.  K(E),  In 
each  component  of  the  transport  medium  are  related  to  the  attachment  and 

stripping  cross  sections  by: 

R,(E)«  -  tu  V  o,(E)«  (10.5] 


“J.R.  Letaw,  J.H.  Adams,  R.  Sllbetberg,  and  C.H.  Tsao,  ■Electron  Capture 
Decay  of  Cosmic  Rays,”  Ap.  Space  Sci.  114,  365*379  (1985) 
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and 


-  n*  V  a.(E)*  [10.6] 

The  electron-captuire  decay  half  lives,  (r^),  and  electron  stripping  (a.)  and 
attachment  (a+)  cross  sections  are  discussed  by  Letaw,  Silberberg,  and  Tsao.“ 
The  electron- capture  decay  halflife  is  approximately  double  the  laboratory 
measurement  because  there  Is  generally  only  one  electron  in  the  K- she  11. 


[10.1]  and  [10.2]  may  be  combined  to  give  an  expression  for  the  total  change  in 
a  cosmic-ray  species  as  a  result  of  electron- capture  decay: 


a[N“,(E)  +  N\(E)] 

at 


R^(E),  N‘,(E) 


R^(E)^  N‘i(E) 


[10.7] 


The  analogous  equation  for  the  cosmic-ray  flux  spectnim  with  pathlength  as  an 
independent  variable  is: 

aiJ^E)  +  J*,(E)] 

-  -  .  r,(E), 


where  the  inverse  mean  free  pathlengths,  r^^,  have  the  same  relation  to  the  rates, 
as  in  [9.5]  and  [9.6].  It  should  be  noted  that  the  summation  in  [10.8]  will 
generally  have  a  contribution  from  only  one  nuclide,  l.e,,  the  species,  j,  which 
can  decay  into  species  1  by  electron  capture. 

In  most  cases  cosmic -rays  achieve  an  equilibrium  charge  state  after  passage 
through  an  amount  of  material  which  is  small  compared  to  other  propagation 
parameters  such  as  the  interaction  mean  free  path  and  the  escape  mean  free  path. 
The  fraction  of  particles  with  a  single  electron  attached  is: 


■I 


J‘,(E)  +)  rJB)^J‘,(E) 


[10.8] 


((r.  -  r,  +  r^)*  +  A  r.  r.]'^  •  (r.  -  r,  r^) 

2  r. 


[10.9] 


“J.R.  Letaw,  R.  Silberberg,  and  C.H.  Tsao,  "Propagation  of  Heavy  Cosmic  Ray 
Nuclei."  Ap.  J.  Suppl.  56,  369-391  (1984) 
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For  each  apecies  a  function  g,(E)  vdiiich  relates  the  flux  of  particles  with  one 
electron  attached,  J‘,  to  the  total  particle  flux,  J®  +  J‘,  at  equilibrium  may  be 
defined  based  on  [10.9]. 


g 


L  J®  +  J‘  Jeq 


[10.10] 


This  function  must  be  obtained  by  proper  weighting  of  the  stripping  and 
attachment  rates  over  the  components  of  the  ISH. 

When  the  attachment  rate  Is  much  smaller  than  the  stripping  rate  plus  the 
electron* capture  decay  rate,  that  Is,  when  only  a  small  fraction  of  cosmic  rays 
have  an  attached  electron. 


g.(E)  -  (S  r*(E)i*)  /  (r^(E),  +  S  r.(E)g,) 

K  X 


[10.11] 


where  the  inverse  mean  free  pathlength  for  attachment,  r+(E),  and  for  stripping, 
r.(E) ,  have  the  same  relations  to  the  cross  sections  as  [8.5]  and  [8.6],  [10.11] 
is  generally  applicable  for  cosmic -ray  transport,  the  exception  being  low-energy, 
ultraheavy  cosmic  rays  in  the  actinide  group. 

Using  the  function  g(E)  to  eliminate  J’,  [10.8]  nay  be  rewritten: 


aJ,<E) 

dx 


’  r„(E).  &<E)  J,(E) 


r«(E)(,  g<(B)  J,(E) 


[10.12] 


Yiou  and  Ralsbeck**  have  demonstrated  that  the  electron  capture  decay  of  ’Be  is 
Important  below  20  MeV  per  nucleon.  Detection  within  the  heliosphere  may  be 
impossible  at  that  energy.  Similar  conclusions  apply  to  secondary  nuclei  such 
as  ”Ar,  "Ca,  *V,  ’’Cr,  and  “Fe  in  the  absence  of  distributed  acceleration 
processes.  Electron  capture  decay  generally  occurs  for  ultraheavy  nuclei. 


“F.  Yiou  and  G.H.  Ralsbeck,  Ap.  J,  7,  U29  (1970) 
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11.  The  Zqueclcm  o£  Distributed  Aecelersticra 

The  transport  equation  described  in  Sections  4  •  10  Is  suoioarized  here.  For 
purposes  of  simplification  we  assume  that  (1)  nuclear  interaction  fragments 
maintain  the  velocity  of  tlieir  progenitors,  the  "strai^t-ahead”  approximation, 
and  (2)  the  interstellar  medium  consists  of  one  substance,  namely  hydrogen. 

The  first  assumption  is  generally  applied  for  studies  of  cosmic  rays  heavier  than 
helium  and  allows 


«0 

I  r,(E,£)^  J,(f)  d£  — >  r,(E)^  J,{E)  [11.1] 

0 

where  r|(E)y  is  related  to  the  partial  cross  sections  and  is  different  from  rf(E)i 
\Aiich  is  related  to  the  total  inelastic  cross  section. 

The  second  assumption  is  generally  applied  because  of  the  lack  of  available 
cross  section  data  for  other  target  nuclei.  Some  investigations  of  the 
importance  of  helium  as  a  target  in  the  interstellar  medium  have  been  performed.” 

The  complete  equation  of  distributed  acceleration  is: 


8J,(E)  Q,(E)  d  r 

-  «  -  .  r.(B),  J,(E)  —  «(E),  J,(E) 

dx.fi  3E  L 

E 

-  r,(E),  J,(E)  4-  (y-Dv'^p-^l  r.(E),  p***  J.CE)  d£ 

0 

-  r,(E),  J,(E)  ^  r,(E)^  J,(S) 

J 


“For  example,  P.  Ferrando,  P.  Goret,  and  A.  Soutoul.  “On  the  Importance  of 
Interstellar  Helium  for  the  Propagation  of  Heavy  Cosmic  Rays,"  ?roc.  19th 
Intemat.  Cosmic  Ray  Conf .  3.  61  (1985) 
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12.  Extausioa  to  an  InhoBOganaoua  Interstallar  Hadlua 

Electron  capture  decay  processes  are  sensitive  to  density  variations  in  the 
Interstellar  nediua.**  Most  o£  these  nuclides  have  a  very  short  aean  lifetime 
compared  to  the  time  during  which  cosmic  rays  move  from  one  density  region  to 
another,  hence  no  effect  on  composition  is  observable.  A  few,  such  as  **Ti,  *‘Nb, 
•^r,  and  *®Ib  are  quite  sensitive  to  possible  density  irregularities.  Measured 
abundances  of  these  isotopes,  particularly  at  low  energy,  may  eventually  offer 
the  opportunity  to  study  aspects  of  the  density  distribution  in  the  ISH  by 
utilizing  cosmic-ray  spectra. 

The  procedure  for  treating  one  type  of  density  inhomogeneity  is  discussed  here. 
Cosmic  rays  are  assumed  to  encounter  density  fluctuations  in  the  ISM  which 
randomly  distributed  in  time  and  space.  These  fluctuations  are  described  by  a 
distribution  function,  F(n).  of  interstellar  gas  number  densities,  n.  The 
function  is  normalized  so  that: 


«0 

I  F(n)  dn  -  1  {12. IJ 

0 

K  rough  estimate  of  this  distribution  function  may  be  made  based  on  models  of 
the  Interstellar  medium,'^  if  the  entire  ISN  is  accessible  to  cosmic  rays. 

The  model  may  be  implemented  by  writing 

dx(n)  «  f(n)  dn  [12.2] 

where  f(n)  is  the  fraction  of  total  pathlength  over  which  the  interstellar  gas 
density  is  between  n  and  n+dn,  and  x<n)  Is  the  portion  of  the  total  pathlength 
with  density  less  than  or  equal  to  n.  It  follows  that 

£(tt)  -  n  F(n)  /  <n>  [12.3] 

where  <n>  is  the  mean  density  of  the  ISK  relative  to  the  distribution  F(a). 
f<n)  will  necessarily  be  normalized  to  unity  If  constructed  in  this  manner. 

For  processes  where  dJ/dx  Is  independent  of  the  density  of  the  Interstellar 
medium,  such  as  nuclear  fragmentation  C6.A],  we  have: 


“J.R.  Letaw,  R.  Sllberberg,  and  C.H.  Tsao,  ■Electron  Capture  Decay  of  Cosmic 
Rays:  A  Model  of  the  Inhomogeneous  Interstellar  Kedlum."  Proc.  19th  Xnternat. 
Cosmic  Ray  Conf.  3.  33  (1985) 

”For  example,  C.F.  McKee  akA  J.?.  Ostriker,  Ap.  J.  218,  I08  (1977) 
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dJ  a  J  dx  <X  J  f(n)  dn  (12.4] 

and  consequently 

la  J  a  1  [12.5) 

because  f(n)  Is  noraallzsd  to  unity.  The  flux  of  cosmic  rays  subject  to  this 
type  of  process  is  Independent  of  the  interstellar  gas  density  distribution. 

For  processes  where  dJ/dx  is  inversely  proportional  to  the  density  of  the 
interstellar  medium,  such  as  nuclear  decay  by  particle  emission  [9.4],  we  have: 

dl  a  J  <n>‘‘  dx  <x  J  n  F(n)  <n>‘*  dn  [12.6] 

and  consequently 

la  J  a  <n>‘  .  [12.7] 

The  flux  of  cosmic  rays  subject  to  this  type  of  process  is  dependent  on  the  mean 
density  of  the  ISH. 

In  the  case  of  electron  capture  decay  ([10.10]  and  [10.11]) 

4J  oc  J  - -  dx  [12.8] 

1  +  kj  <n> 

where  k,  and  k,  are  energy-  and  species -dependent  functions  which  are  independent 
of  the  ISM  density  distribution.  In  the  two  extreme  cases,  k,  «  1  and  »  1, 
the  solution  to  this  equation  is  Identical  In  form  to  [12.5]  and  [12.7], 
respectively.  For  intermediate  cases,  the  solution  Involves  an  integral  over 
the  density  distribution  of  the  Interstellar  medium  which  may  provide  information 
about  inhomogcnelty  over  the  cosmic-ray  path. 
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